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Abstract: Robot motor capability is a crucial factor for a robot, because it affects how 
accurately and rapidly a robot can perform a motion to accomplish a task constrained by 
spatial and temporal conditions. In this paper, we propose and derive a pseudo-index 
of motor performance (pip) to characterize robot motor capability with robot kinematics, 
dynamics and control taken into consideration. The proposed pl^ provides a quantitative 
measure for a robot with revolute joints, which is inspired from an index of performance 
in Fitts's law of human skills. Computer simulations and experiments on a PUMA 560 
industrial robot were conducted to validate the proposed pl^ for performing a motion 
accurately and rapidly. 

Keywords: robot motor capability; pseudo-index of motor performance; speed-accuracy 
constraint; Fitts's law 



1. Introduction 

Robots are widely-used in manufacturing tasks, such as assembly [1], cutting [2], deburring [3], etc. 
Without loss of generality, these manufacturing tasks can be decomposed into a sequence of motor tasks 
that are usually described by spatial and temporal constraints on objects in the environment [4-6]. Thus, 
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a major problem for a robot to accomplish a task is to satisfy the spatial and temporal constraints of a 
task. A spatial-temporal constraint requires a robot to perform motor movements accurately in a desired 
trajectory task. Usually, a robot is asked to achieve a given position with a minimum movement time. 
However, each robot has its own motor capability according to its mechanical and electrical system and 
may not be qualified for tasks with strict spatial-temporal constraints. 

Robot motor capability is the motor ability for a robot to accomplish a task constrained by spatial 
and temporal conditions. Traditionally, a robot was evaluated by its accuracy and repeatability [7-10]. 
American National Standard defined the static and dynamic performance characteristics for industrial 
robots and robot systems [11,12]. Throughout the previous work, there has not been a performance 
metric to measure the robot motor capability of a robot motor system for tasks constrained by spatial and 
temporal conditions. Thus, there is a need to develop a quantitative measure of a robot motor system. 

In this paper, the proposed method for characterizing the motor capability of a robot is inspired by 
Fitts's law, which is one of the well-known metrics for studying human rapid movements. Fitts's law was 
revealed by Paul Fitts in 1954, showing the information capacity of a human motor system. In Fitts's 
law, the motor performance is the ability to consistently produce movements and is described by two 
main factors — speed and accuracy. Fitts obtained the following equation, often called Fitts's law, as: 

2D 

T^, = a + 6.1og2(— ) (1) 

where Tmt is the minimum movement time and a and b are proportional constants, and Fitts quantitatively 
defined an index of task difficulty (Id) as log2(^) (bits/response) to specify the minimum required 
information for rapidly-aimed movements to accomplish a given task, where is a target distance and 
Vl^ is a target width defined as the allowable tolerance along the moving direction. Fitts also defined an 
index of performance, 7^, to indicate the fixed information capacity: 

1 2D 

and found that Ip is relatively constant over a considerable range of movement amplitudes and 
tolerance limits. 

Actually, Fitts's equation is a special case of the quasi-power function with a varying exponent, n, as 
n ^ oc [13]. The quasi-power function with exponent n is shown as T^t = + &n • (^)^- When 
n ^ 1, the quasi-power function becomes a linear function as: 

Tmt = ai + h-{^) (3) 
Equation (3) is re- written as a linear speed-accuracy relationship: 

W = a+b'-{^) (4) 

-l-mt 

where a and are proportional constants, and the index of performance (Ip) of Fitts's law in Equation (3) 
is defined as ^ by ignoring the constant, ai. Thus, ^ in Equation (4) indicates the index of performance, 
as well, because smaller bi results in smaller T^t and smaller b' if D and W are fixed in both equations. 

The paper proposes the method to derive the speed-accuracy constraint of a robot system the same as 
Equation (4), and therefore, the pseudo-index of motor performance (pip) of the robot system is defined 
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as ^ in Equation (4). The speed- accuracy constraint is characterized with robot kinematics, dynamics 
and control taken into consideration. Thus, the mathematical derivations of pip are applicable to other 
robots with revolute joints. Furthermore, pip is a quantitative measure for the motor capability of a 
robot. The basic concept of our approach is to treat a robot motor system as an information channel 
and determine the robot motor capability. Since the information capability of a channel is affected by 
the inaccuracy of kinematics and dynamics models, we can quantitatively measure it by the variability 
of movements that a robot aims to produce. To verify the proposed method, computer simulations and 
experiments on PUMA 560 and 260 industrial robots are conducted to verify the validity of the proposed 
approach in understanding the motor capability of a robot motor system. 

This paper is organized as follows. In Section II, we formulate a robot task as an optimization 
problem constrained by the proposed speed-accuracy constraint derived from the kinematics, dynamics 
and control of a robot motor system. In Section III, we present the proposed pseudo-index of motor 
performance (pip) derived from the speed- accuracy constraint and show how pip is utilized to evaluate 
the motor capability of a robot. In Section IV, computer simulations and experiments on PUMA 560 and 
260 industrial robots are presented. Discussions and conclusions are summarized in Section V. 

2. Speed- Accuracy Constrained Optimization 

2.7. Problem Formulation 

We formulate a robot task as an optimization problem subject to spatial and temporal constraints. 
The proposed speed-accuracy constraint characterizes the relationship between the robot joint speed and 
the Cartesian position error of the end-effector of the robot. By utilizing the proposed speed-accuracy 
constraint, the optimization problem is solved without violating the robot motor capability. We first 
formulate a basic rapid robot movement that moves its end-effector to reach a desired joint location, 
0 = [^1, ^2, • • • i^ii' ' ' 1 ^nY ^ with an allowable tolerance, e^, Cy and (i.e., spatial constraint), in 
the Cartesian space by the maximal linear speed (i.e., temporal constraint), where is the number 
of degrees of freedom (DOF) and the superscript, T, denotes a matrix transpose. Here, we need to 
determine the joint velocities that introduce the maximal linear speed of the end-effector. The joint 
velocities are represented by a vector as ^ = [^i, ^2, • ' ' ^^ir ' ' , ^a^]^. and the basic rapid movement of 
the end-effector of a robot is formulated as an optimization problem as: 

max||J^(6l)^||2 

subject to \dx\ < e^, \dy\ < Cy, \dz\ < e^, and 

0 < 1^,1 < forz = 1 to (5) 

where || • II2 is the Euclidean-norm operation, J^{0) is the 3 x Jacobian matrix relating the joint 
velocity to the Cartesian velocity of the end-effector, dx, dy and d^ are the respective Cartesian position 
errors, e^,, Cy and > 0 are the allowable position errors and 9^^^^ > 0 is the maximum joint, velocity 
that can be driven. Obviously, 0^^^ in Equation (5) may not be the optimal solution to the problem, 
because it may violate the spatial constraint described by Cx, Cy and e^. Thus, both the spatial and 
temporal constraints must be satisfied simultaneously in this optimization. To consider these constraints, 
an appropriate speed-accuracy constraint must be derived. 
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2.2. Robot Speed-Accuracy Constraint 

Since the spatial inaccuracy of the end-effector of a robot with revolute joints is caused by the 
inaccuracy of robot kinematics and dynamics models and the disturbance of the environment [8], we 
represent the Cartesian position error of the end-effector of an A^-DOF manipulator (e.g., a PUMA 
robot) with revolute joints by a linear model expressed as: 



4 



= [J^{e)Ad + W Aa + J- W Aa + jf (0) A^] + J^{e)Ae 



c{e) + JM^e 



(6) 



where dx, dy and d^ are the respective Cartesian position errors along each coordinate axis, J^{0), 
J^{0), J^{6) and J^{0) are the 3 x kinematic-error matrices [8], J ^{6) is the 3 x Jacobian 
matrix, 0 is an x 1 joint position vector, (d, a, a, (3) are the D-Hparameters [14] in an A^ x 1 vector 
form, (Ad, Aa, Aa, A/3) denote small changes in the corresponding parameters and C{0) is a 3 x 1 
vector, summing all the kinematic errors, due to the D-H parameters (d,a, a,/3); this is due to the 
mechanical tolerance in manufacturing the links and joints of the robot. 

To derive /\0 in Equation (6), we assume that the widely-used "computed torque" technique with 
a proportional-plus-derivative (PD) controller [15] is used. Thus, each robot joint is modeled as a 
single-input-single-output system with disturbance from other joints. The coupling effects from other 
links are considered as a disturbance to the robot. A PD control is utilized to servo the manipulator 
system to track the desired trajectory, and the feedforward torques along the actual trajectory are 
computed to compensate for the nonlinear disturbance. A block diagram showing the PD control scheme 
for joint, z, is shown in Figure 1. 

Figure 1. Proportional-plus-derivative (PD) control scheme of a robot joint, i. 




lis) 



In Figure 1, Jij{0), /f, G^{0) and Vf are, respectively, the computed counterparts of the actual 
effective inertia terms, Jij{9), effective damping coefficient, gravity terms, G-{0), and Coulomb 
friction, V^, at the actual shaft of joint, i, of the manipulator, is the gear ratio, armature resistance, 
motor constant, K^, and the back electromagnetic force constant, Kl, is known dc motor parameters 
for joint, i. K'l and Kl are, respectively, position and velocity feedback gains of the controller. Their 
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values are chosen subject to the mechanical resonant frequency constraint of the manipulator. Of{s) and 
9i{s) are, respectively, the desired and actual angular displacement of joint, i. T}^{s) is the disturbance 
in Laplace transform and T^jf{s) is the computed feedforward torque compensating for the disturbance. 

In order not to excite the mechanical resonant frequency of the manipulator, the undamped natural 
frequency is set to no more than one-half of the structure resonant frequency. Complying to this 
constraint, we can obtain the following relation and the upper bound of Kp. 

0<K'< ^Q^'^Q^'^^ (7) 

a 

where uUq and Jq are, respectively, the measured structural resonant frequency and inertia of joint, i, at a 
known location. Since Figure 1 is a second-order system, the damping coefficient, Q, of the system can 
be obtained as: 

^ B}J, + KIKI + KIKI 
In order to avoid the oscillatory underdamped response, Q is set to >1. Thus, becomes: 

V — jy-i W 

a 

Because the actual Ju and fi are usually not available, we use their computed counterparts, Jf- and , 
from the manipulator dynamic model. Substituting Equation (7) into Equation (10), we obtain Kl as: 



= (10) 



To determine AO in Equation (6), we utilize Figure 1 to evaluate the steady-state error of 0 by 
a velocity-related joint position command, 9d{t) = At, where A is the amplitude of the input. The 
steady- state error of 9i is derived as: 



K'K 

a p 



e AHi{d)e + AGi{d) + AVi 



(11) 



where AH^^,{e) = H^^,{e) - Hf^,{0), AG,{e) = 0,(0) - Gl{d\ AV, = V,- and AH^ is an 
N X N symmetric matrix, whose elements are Ai/-^^(0), where j and k are integers and 1 < j < 
and 1 < A: < iV. 

Replacing the elements of Ad in Equation (6) with AOi in Equation (11), we obtain the Cartesian 
position errors of the end-effector of the robot in the form of speed-accuracy constraint as: 



dx 

dy 

d,. 



T 



c{d) + j.{d)KAP{d) + JJO)K^ AH{d)e 



(12) 
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where AP{0) is an x 1 vector and AP{e) 



AG(6>) + AV, where AG = AG ^{6) 



(AG.ie) + AVi, ■ ■ • , AGj,{9) + AVn 

1 T 



^AG^id) 



N X N diagonal matrix, whose elements are 



4n^ 



and AV = 

4n^ 



[AV^i,--- ,Al^^f,li:isan 



AH{d) is an A^^ 



X 



matrix, and AH{e) = diag{Affi(0), • • • AffN(^)}, ^ is an x N matrix and 

, 0; • • • ; 0, • ■ ■ ,0, 9j^], where = ^ for i = 1 to N, and 0 is an 



, e'^f, where 0- =0^ fori = l to A^. 



* = [0,,O,- - ,O;O,02,O, 
A^^ X 1 matrix and 0 = [^^ , 

From Equation (12), it shows that when the joint speed increases (Onew = i^^oid, ^ > 1) and the 
desired joint positions, 6, are fixed, a larger absolute Cartesian position error (\dx\, \dy \ and \dz\) will be 
generated. Thus, Equation (12) indicates the trade-off relationship between the spatial accuracy and the 
joint speed. 



3. Pseudo-Index of Performance for a Robot System 

The purpose of determining the pseudo-index of the performance, pip, of a robot motor 
system is to obtain a quantitative measure. After the speed-accuracy constraint is obtained, 
we can derive and calculate the index of the performance, pip, of a robot motor system 
from this constraint. Thus, we define the pseudo-index of the performance (pip) of the linear 
speed-accuracy equation in Equation (4) as pip = ^ to quantitatively indicate the amount of 
information capacity of a robot motor system. According to the speed-accuracy constraint (see 
Equation (12)), C{6) + J^{6)KAP{6) can be referred to as a' in Equation (4). The joint velocity 
command is expressed as 6d = pOo, where is a reference joint velocity vector and p > 0 
is a scale factor to change the joint velocities. We can express Equation (12) in the form of 
Equation (4) as: 



dx 

dy 



[c{e) + j^{e)KAP{e)] + 



1 

pipx 

1 

pipy 
1 

L pipz 



' P 



(13) 



where pipx, pipy and pIpz are the pseudo-indices of performance on the x-, y- and z-axis of the Cartesian 
space, respectively, and are derived as: 

-1 



pi,, = [IJJl0)K%AH{0)@o 



pipy 
PIpz 



Iy'J^{0)K%AH{0)@o 
I/J^{0)K^IaH{0)@o 



(14) 



where Ix, ly and Iz are the unit vectors along the x-, y- and z-axis, respectively, #o is an A^^ x matrix 
and #0 = [^1, 0, • • • , 0; 0, 0> ' ' » 0; ' ' ' ! 0) • ■ ,0, 0j^], where 0^ = for « = 1 to A^, and 0o is an 
xl matrix and ©o = [^i , • • • , Oj^f, where 0- = for i = 1 to A^. 
From Equation (13), we obtain the quadratic function of the Cartesian position errors, (dx, dy, d^), 
with respect to the scalar of the joint velocity, p. Figure 2 illustrates the obtained quadratic function of 
the Cartesian position error of the end-effector of a robot along the x-axis (dx) with respect to the scalar 
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of the joint speed (p), where the intercept on the axis of the Cartesian position error of the end-effector 
of a robot is Ix^[C{6) + J^{6)KAP{6)]. From Figure 2 and Equation (13), it is shown that a robot is 
capable of performing more accurately with a larger p/^ when the joint velocities, Od, are fixed (i.e., p is 
fixed). In other words, when the Cartesian position tolerance, e, of a given task is fixed (i.e., is fixed), 
a robot with a larger pl^ is capable of accomplishing the task with faster joint velocities. Interestingly, 
the pseudo-index of performance is advantageous to show the trade-off relationship between speed and 
accuracy, if we assume that there exists a fixed amount of information capacity of a robot motor system. 

Figure 2. Quadratic function of the Cartesian position error of the end-effector of a robot 
along the x-axis (d^^) with respect to the scalar of the joint velocity (p), where the quadratic 
coefficient is {plpx)~^ and Od = pOo. 

d (Cartesian Position Error Along the X-axis) 



' 

(Scalar of Joint Speed) 

4. Computer Simulations & Experimental Work 

Computer simulations were performed by MATLAB Toolbox [16], and experiments were performed 
on 6-degree-of-freedom PUMA 560 and 260 robots [14] with revolute joints to validate the proposed 
speed-accuracy optimization and demonstrate that the motor capability can be measured by the proposed 
pseudo-index of performance (pip) on a rapid-moving task. In the experiments, the Cartesian position 
error was with respect to the 2D space, because we followed Fitt's experiments, where the position error 
was only measured in the 2D space. The Cartesian position error on the last dimension (z-axis) would 
not affect that on the other two dimensions (x- and y-axis), because they are independent. 

4.1. Simulation on a PUMA 560 Robot 

Figure 3 shows the scheme of our task. The robot was asked to move its end-effector to hit a target 
position (G) and stop at an end position (E) by various configurations of movement speed (D/T^t) and 
distance (D). With each configuration of movement speed and distance, this resulted in a Cartesian 
position error (Vl^) at the end-effector. 

In this simulation, a PUMA 560 industrial robot was asked to move with the maximum speed to hit the 
target position within the spatial tolerances, ax and ay, which denote the maximum position tolerances 
along the x- and ^-axis, respectively. For simplicity, the robot only used joints 1, 2, and 3 to perform 




r[C(0) + JA0)KAP(0)] 
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the task. Based on the task description, we formulated the task as an optimization problem subject to the 
spatial-accuracy constraints for maximizing the robot speed: 

max II J^(^)^||2 

subject to 141 < ^x, \dy\ < (Jy, and 0 < |^^| < ^f^^ (15) 

where JjyO) is the Jacobian matrix, 4 and dy are the Cartesian position errors of the end-effector of the 
robot along the x— and t/— axis, and 6^^^ is the maximal velocity of joint, i {i = 1,2, 3), that can be 
driven. The parameters of the PUMA 560 robot were from [17]. We set the distance from the tip of the 
pin to the coordinate frame of the sixth joint to 0.1 m. 

Figure 3. Scheme of our task, where the moving position error is defined as the error between 
the actual and target positions. (S: start position, G: target position and E: end position). 




End position ^ 



First, we evaluated how the Cartesian position error was affected by Am, AG{0) and AV{0) under 
a fixed joint distance. We used joint-interpolated motion planning for the robot. Simply, we subtracted 
C{6) from the Cartesian position errors of the end-effector of the robot. The joint distance, Odiff, 
between the initial and target joint positions was [7r/4, 7r/4, — 7r/4, 0, 0, 0], and the execution time for 
all the joints was set to T^^. Thus, the joint velocities were represented as Odiff/Tmt (rad/sec). If we 
set ^0 as the nominal joint velocities finishing Odiff by one second, the joint velocities were expressed 
as p • where p = l/|T^t|. The joint velocity, acceleration and the inverse dynamics of the robot 
were calculated by the recursive Newton-Euler equations [14,18], and the Cartesian position error of the 
end-effector of the robot was calculated by J^{0)AO and verified by Equation (12). 

When the errors between the computed and actual masses. Am ^rri^ ^ 

(±0.1%), (±0.2%), • • • and (±1%), were assumed as the cause of the Cartesian position errors, we 
simulated four cases of a combination of A V (Coulomb friction) and AG{0) (errors between the actual 
and computed gravity terms). Their results are explained and discussed as follows: 

Case (I) AV = 0 and AG {6) = 0: In this case, the simulation results showed that the Cartesian 
position errors of the end-effector of the robot (dx = dy = d^ = 0) almost do not exist. Apparently, they 
can be validated by Equation (12), because of AP{e) = 0 (AV = 0, AG{e) = 0 and AH{e) = 0). 

Case (2) AV ^ 0 and AG {6) = 0: In this case, the simulation results of the Cartesian position 
errors along the x- and ^-axis (dx and dy) were shown by the blue dots with Am = 0 in Figure 4. The 
Cartesian position errors were only due to A V of Equation (12). 

Case (3) AV = 0 and AG{0) ^ 0: In this case, the simulation results showed that when p = 1/2 
and Am was changed from 1.0% to —1.0%, the Cartesian position errors were shown by the line of 



Sensors 2013, 13 



8420 



red crosses in Figure 4 and caused by AG{0) and AH{0) of Equation (12). On the other hand, when 
the joint velocities were changed by p (1/2, 1/4, 1/6), Figure 5 shows that the slope of the line showing 
the Cartesian position errors was also changed. It shows that as p= 1/6 (gray line, lowest speed), the 
Cartesian position errors were least (highest accuracy). In addition, when Am was fixed. Figure 5 
shows that the Cartesian position errors caused by different p were in a linear line (Gi pink line), and 
the increment of the Cartesian position errors was proportional to p^, because J^{0)K^ AH{9)& in 
Equation (12) indicates that the Cartesian position errors were proportional to {0 ) (i.e., p^) when Am 
was fixed. However, none of the pink fines in Figure 5 did passed through the origin, but a point G^, 
(i = • • • , —2, —1, 1, 2, • • • ). These results were clearly shown by Equation (12), because the Cartesian 
position errors resulted from AG{0) and AH{6), but when the joint velocities were set to zero, the 
Cartesian position errors were only due to AG{0) {i.e., Gi). 

Case (4) AV ^ 0 and AG{0) ^ 0: This case has the same results as Case 3, except the lines 
showing the Cartesian position errors caused by fixed p did not pass through the origin. 

Figure 4. Simulation results of the Cartesian position errors of the end-effector of a PUMA 
560 robot, (a) Case 1: AV = 0 and AG{e) = 0 (red dot); (b) Case 2: AV ^ 0 and 
AG{e) = 0 (blue dot); (c) Case 3: AV = 0, AG{e) 7^ 0 and p = 0 (red crosses); 
(d) Case 4: AV 7^ 0, AG{e) 7^ 0 and p = 0 (blue crosses). 
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Then, we evaluated the pip of a PUMA 560 robot with various Am's and Cartesian distances, where 
Am = m — m^, m were the actual masses and m^ were the computed masses. In other words, 
a robot with a different Am had a different motor capabifity, because Equation (14) shows that the 
pseudo-index of performance is affected by AH [6) (i.e., Am). We set the Cartesian distances between 
the start (5^) and the target (G) positions to 0.0508, 0.1016, 0.2032 and 0.4064 m. The four movement 
velocities (D/Tj^t m/sec) for a specified distance (D) were set in Table 1. To investigate the Cartesian 
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position errors of the end-effector of the PUMA 560 robot caused by the changes of the joint velocities, 
we subtracted the velocity-independent terms, C{6) + J^{6)KAP{6), from the Cartesian position 
errors at the target position. For D = 0.0508, 0.1016, 0.2032 and 0.4064 (m), the joint velocities were 
represented as (0^jj)/T^t= p • where is the nominal joint velocities finishing Odiff by one 
second, p is a scalar and its value is p = l/|r^t|. We changed the joint velocities by changing p with 
a 0.01 incremental value. The joint velocity, joint acceleration and the inverse dynamics of the robot 
were calculated by the recursive Newton-Euler equations [14,18], and the Cartesian position error of the 
end-effector of the robot was calculated by J ^{6) AO. 

Figure 5. Simulation results of the Cartesian position errors of the end-effector of a PUMA 
560 robot in Case 3 when the joint velocities are non-zero: AV = 0, AG{6) ^ 0 and 
p ^ 0 (crosses). 
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Table 1. Four movement velocities, D/Tmt m/sec under different D' s. 

D (m) D/T^t (m/sec) 

0.0508 0.0508/0.392 0.0508/0.281 0.0508/0.212 0.0508/0.180 

0.1016 0.1016/0.484 0.1016/0.372 0.1016/0.260 0.1016/0.203 

0.2032 0.2032/0.580 0.2032/0.469 0.2032/0.357 0.2032/0.279 

0.4064 0.4064/0.731 0.4064/0.595 0.4064/0.481 0.4064/0.388 



Figure 6 shows the Cartesian position errors along the x-axis {d^} with respect to D /Tmt when 
D = 0.0508 and Am was (a)— 1%, (b)— 0.6% and (c)— 0.2%, respectively. In each sub-figure, we 
used a quadratic function to do curve fitting for our simulation results. Since the quadratic coefficient is 
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the reciprocal of pipx, by referring to Figure 2, we noticed that Figure 6(c) had the greatest pipx among 
Figure 6(a-c). It was explained by the reason that Figure 6(c) had the least amount of Am = —0.2%. 

Figure 6. Simulation results of the Cartesian position errors of the end-effector of a PUMA 
560 robot along the x-axis (dx) with respect to the Cartesian velocities (D/Tmt) when 
D = 0.0508 and Am was (a) -1%; (b) -0.6%; (c) -0.2%. Red triangles represent the 
simulation results; the blue curve represents curve fitting. 




5 6 7 8 9 10 11 12 



2.54*D/Tmt (m/sec) 
(c) 

To verify the robot had different motor capabilities with respect to p on the x-axis, we changed Am, 
fixed D (i.e., D = 0.0508,0.1016,0.2032 or 0.4064) and compared the quadratic coefficients among 
these various Am's. Figure 7 shows the quadratic coefficients when Am was —1%, —0.6% and —0.2% 
in each sub-figure, where (a), (b), (c) and (d) represent D = 0.0508, 0.1016, 0.2032 or 0.4064. We found 
that the robot with Am = —0.2% had the largest pIpx compared to Am = — 1% and —0.6%. In 
other words, the robot with the smaller error between actual and computed masses had the better motor 
capability. The simulation results validated Equation (14) because (^ij^ J ^{6)K^^ /\H {6)&^ and 
AH{0) is definitely affected by Am. 
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Figure 7. Simulation results of the Cartesian position errors of the end-effector of a PUMA 
560 robot along the x-axis (dx) with respect to the scalar of the joint velocities (p) when Am 
was —1% (black dot-dashed line), —0.6% (red dashed line) and —0.2% (blue solid line). 
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4.2. Experiments on PUMA 560 and 260 Robots 

To validate the speed-accuracy constraint, we performed the experiment on a PUMA 560 robot (see 
Figure 8). The joint distance, Odiff, between the initial and target positions of a joint-interpolated 
motion planning was [7r/4,7r/4,7r/4, 0,0,0], and the execution time for all the joints was set to 
Tjnt = 1-6, 1.07, 0.8 and 0.64 seconds for tests 1, 2, 3 and 4, respectively. Thus, the joint velocities 
of these tests were 6 = 6diff/Tmt = P ' ^o(rad/sec), where = ^diff/^ are nominal joint velocities 
finishing Odiff by one second and p = l/|T^t| = 0.625, 0.935, 1.25, 1.5625 for tests 1 to 4, respectively. 
The relationships of joint velocities between test 1 and tests 2, 3 and 4 were obtained as 62 = 1.5^i , 
63 = 201 and 64 = 2.5^i. Although we could not obtain the actual masses of the links, we performed 
various Am's (the ratio of the error between the computed masses and the original computed masses we 
had to the original computed masses) by —3%, —2%, —1%, 0%, 1% and 2% and obtained similar results 
in Figure 5. Figure 9 shows the Cartesian position errors (dx.dy) of the end-effector of the PUMA 
560 robot caused by the changes of joint velocities and Am. In Figure 9, each line was plotted by 
the Cartesian position errors of the end-effector of the robot introduced with various computed masses 
and the same joint velocity. This shows that the Cartesian position error, oc Am, is under the same 
joint velocity. 
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Figure 9. Experimental results of the Cartesian position errors of the end-effector caused 
by the changes of joint velocities (p) and computed masses, where Am = —3%; •: 
Am = -2%; ■: Am = -1%; *: Am = 0%; x: Am = 1%; A: Am = 2%. 
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From Figure 9, we obtain the relationship of Cartesian position errors caused by fixed Am and the 
changes of the joint velocities (p), as shown in Figure 10. Figure 10 represents Am for —3%, —2%, 
— 1%, 0%, 1% and 2%, respectively. The ratios of the Cartesian position errors of the end-effector 
of the PUMA 560 robot caused by 62, 63 and 64 to the Cartesian position errors caused by 61 in 
Figure 10 were shown in Table 2. These ratios of errors were close to the expected ratios, 
p2 ((1.5)2 = 2.25, (2)2 = 4, (2.5)2 = 6.25), since O2 = 1.5^i , Os = 20i and 64 = 2.5^i. This 
shows that the Cartesian position error, oc p'^, is under the same Am. Equation (13) explains this 
phenomenon when PUMAS 60 had a fixed plpx- 
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Figure 10. Experimental results of the Cartesian position errors of the end-effector caused by 
the fixed Am and the changes of the joint velocities (p). (a) Am = —3%; (b) Am = —2%; 
(c) Am = -1%; (d) Am = 0%; (e) Am = 1%; (f) Am = 2%. 
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To demonstrate that the pseudo-index of performance (pipx) is useful to evaluate the motor capability 
of various robots with revolute joints, we compared the motor capabilities of the PUMA 560 and 260 
robots by measuring their pIpx, where the PUMA 260 was a small robot with a similar mechanism to 
the PUMA 560. Figures 11 and 12 show the PUMA 560 and 260 robots' quadratic coefficients when 
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D = 0.0508, 0.1016, 0.2032 and 0.4064. As plp^ is the reciprocal of the quadratic coefficient, it is 
obvious that the PUMA 260 robot has better motor capabihty than the PUMA 560 robot. 

Table 2. Ratios of the Cartesian position errors of the end-effector of the PUMA 560 robot 
caused by 62, 63 and 64 to the Cartesian position errors caused by 61. 
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Figure 11. PUMA 560 experimental results of the Cartesian position errors of the 
end-effector along the x-axis {dx) with respect to L^/T^^ (a) D = 0.0508 m; (b) D = 0.1016 m; 
(c) D = 0.2032 m; (d) D = 0.4064 m. 
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Figure 12. PUMA 260 experimental results of the Cartesian position errors of the 
end-effector along the x-axis (dx) with respect to L^/T^^ (a) D = 0.0508 m; (b) D = 0.1016 m; 
(c) D = 0.2032 m; (d) D = 0.4064 m. 
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5. Discussions and Conclusions 

The proposed speed-accuracy constraints showed the relationship between Cartesian position errors 
of the end-effector of a robot and joint velocities. This relationship can be illustrated by a vector 
representation, as shown in Figure 13. For simplicity, Figure 13 shows the 2D projection of a 3D 
Cartesian position error caused by joint position errors (AO). The Cartesian position error on the 
2D plane, AEaii, is composed of dx along the x-axis and dy along the y-axis. It is shown as the 
summation of AEy, A Eg and AEs, where AEy represents the Cartesian position errors introduced by 
the errors between the actual and computed Coulomb frictions, AEq represents the Cartesian position 
errors introduced by the errors between the actual and computed gravity terms and AEs represents 
the Cartesian position errors regulated by the velocity changes. Here, AEq is affected by Am and 
AEg oc where Am = amTn and is a scalar. Also, AEs is affected by 0 and AEs • (0^5)^, 
where ^ = o^s^o. ^0 is a reference speed and ag is a scalar. The angle, 0, is determined by 0, m and 
Am. When the desired joint positions, 0, and link masses, m and Am, are determined, we can choose 
as (maximum speed) based on the vector relationship to make the vector, AEaiu stay in the boundary of 
the task constraints. 
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Figure 13. Vector representation of Cartesian position errors. 




From Equation (14), it is obviously shown that the pseudo-index of performance is the same with 
respect to the joint velocities, 6, instead of the Cartesian speed, dD/dt. Thus, the pseudo-index of 
performance was not the same for the sub-figures of Figures 1 1 and 12. In addition, when a robot moves 
to a target in a fixed distance. Am of the robot will influence its motor capability. The less error in 
masses will give the robot more information capacity for trading off movement speed and accuracy to 
accomplish a task. 

The pseudo-index of performance is also affected by the target joint positions, 0, because from 
Equation (14), J^{0) dramatically affects the Cartesian position errors, dx, dy and d^. Figure 14 shows 
the L^-norm of I^^ • J^{0), ly^ • J^{0) and • J^{0) with respect to index number, where I^, 
ly and Iz are unit vectors along the x-, y- and z-axis, respectively. L^-norm indicates the maximum 
Cartesian position errors, dx, dy and rf^, to which J^{0) can contribute. The index number is generated 
from the data format as index (^jiyAndex(^j2)\index(^j'^yAndex(^j4^^^^^ where index (^ji) is a 

sampling index by joint resolution within the joint range of joint, i [14], joint resolution is determined 
arbitrarily, index (^ji) are the most significant bits (MSB) and index i^jq) are the least significant bits (LSB). 
The index number represents the joint positions for all joints. To illustrate the concept and simphfy the 
computation, joints 4, 5 and 6 are fixed as zero-degree, and the index number is only generated from 
joints 1, 2 and 3 with 20-degree joint resolution. Figure 14 shows the distribution of the L^-norm on 
dx, dy and dz with respect to the index number (joint positions). For some cases, an index number may 
have entirely different degrees of Cartesian position errors along the x-, y- and z-axis. For example, the 
index number, 1,639 (joint 1: 0 degree; joint 2: —125 degrees; joint 3: —45 degrees), has the maximum 
L^-norm on dx, but a relatively small L^-norm on dy. However, this index number is on the range 
boundary of joint 3, and it apparently is not an operation joint position. Similarly, another L^-norm on 
dx whose index number is 1,556 (joint 1: —20 degrees; joint 2: 35 degrees; joint 3: —25 degrees) is not 
on the joint range boundaries. For this index number, dx is also relatively larger than dy. 

In this paper, we have developed and presented a quantitative measure, the pseudo-index of 
performance (pip), to characterize the motor capability of a robot and the speed-accuracy constraints 
to optimize its motor performance for performing a given task. By using the pip, one clearly understands 
the motor limitation of a robot. Since the pseudo-index of performance is based on the dynamics 
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and kinematics of a robot system with revolute joints, the derived equations of the pseudo-index of 
performance are apphcable to other robots with revolute joints. The PUMA 560 and 260 were revolute 
robots and feasibly the examples to validate the pseudo-index of performance. Computer simulations 
and experiments on PUMA 560 and 260 industrial robots have validated the characteristics of the 
speed-accuracy constraint and the pip and how the pip is utilized to measure the motor capability of 
a robot system. 

Figure 14. Cartesian position errors with respect to index number. Joint resolution is 
20 degrees, (a) dx', (b) rf^; (c) rf^; (d) the same index number (joint positions) results in 
different L^-norm values on and d... 
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